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ABSTRACT 

The definition of homology for single- component galaxies is clear, but for multi- 
component (luminous and dark matter) galaxies there is some ambiguity. We attempt 
to clarify the situation by carefully separating the different concepts of homology that 
have been used to date. We argue that the most useful definition is that a set of galax- 
ies is homologous if they are the same in all respects up to a set of three dimensional 
scaling constants which may differ from one galaxy to the next. Noting that we are 
free to choose the dimensional constants, we find that a set of hydrodynamic simulated 
galaxy merger remnants is significantly closer to homologous when the dimensional 
length constant is taken to be the radius containing equal amounts of dark and bary- 
onic matter rather than the usual observationally motivated choice of the baryonic 
half-mass radius. Once the correct dimensional scaling constants are used, the stel- 
lar velocity dispersion anisotropy is essentially the sole source of the variation in the 
kinematic structure of these simulated merger remnants. In order to facilitate the use 
of these scaling constants to analyse observed galaxies, we calculated the relationship 
between our preferred dimensional scaling constants and the typical observationally 
accessible quantities. 



1 INTRODUCTION 



To what extent are galaxies (or some subset of galaxies) ho- 
mologous in the sense that all individual objects are scaled 
copies of one another? The discovery of the Fundamental 



Plane (FP) ( |Djorgovski fc Davis|1987||Dressler et al.|1987[ ) of 
elliptical galaxies distinctly sharpened this question because 
galaxies showed a regularity that almost, but not quite, fol- 
lowed that expected from the virial theorem for systems in 
dynamical equilibrium. The virial theorem holds regardless 
of whether galaxies are homologous, but if each galaxy were 
unique, then the FP would have a large scatter. If all galax- 
ies were simply scaled copies of one another, then the FP 
would follow the virial expectation. In either of these cases, 
the FP would likely not have generated the interest that 
it has because it would either be less useful (in the former 
case) or less interesting (in the latter case). Actual galaxies 
seem to be in between these two possibilities: they are quite 
regular in structure, but do not quite follow the simple virial 
expectation. An understanding of the origin of the FP re- 
quires at least one ingredient aside from the virial theorem, 
and the search for that missing ingredient has been the focus 
of a great number of observational and theoretical studies. 

Homology means that all galaxies are somehow scaled 
copies of one another. Obviously the property of homology 



applies only to a set of objects, not to individual galaxies. 
If galaxies were made up of only one type of matter, then 
there is no ambiguity in the definition of homology. However, 
they are in fact made up of dark matter, luminous matter, 
and gas, the behaviour of each component is dominated by 
different physics, and the components are only weakly cou- 
pled. For such systems, the proper definition of homology is 
not clear. Several informal definitions seem to be in common 
use; for single component systems, all of these definitions are 
equivalent. Not so for multi-component galaxies. 

We formalize the definitions currently in use and dis- 
cuss which definitions are stronger or weaker in the sense 
that the galaxies that fit one definition are or are not a 
subset of those that fit another definition. Then we apply 
these ideas to galaxy remnants produced in numerical sim- 
ulations. Finally, we discuss how these ideas may be applied 
to observed galaxies. 

It is obviously the case that homology can only ap- 
ply to some subset of galaxies (ellipticals and spirals are 
quite distinct). It has also been established that even for re- 
stricted sets of galaxies, homology can only be approximate 
because, for example, the shape of the surface brightness 
as a function of radius (parametrized by the Sersic parame- 
ter n) correlates with galaxy luminosity ( Prugniel & Simien 
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1997). However, elliptical galaxies are still quite close to be- 



ing homologous-close enough that the concept is still very 
useful. 

In Section [2] we show that homology is closely related 
to the choice of units for each galaxy: homologous galaxies 
are all the same up to a choice of units, where the units are 
allowed to change from one galaxy to the next. For multi- 
component galaxies, there is an ambiguity in the definition 
of homology that we may be able to exploit: perhaps the 
baryonic half-light radius is not the right choice to make all 
galaxies look the same. In Section [3] we review recent ob- 
servational and theoretical work on the subject of galaxy 
structure and homology. In Section [4] we use numerical sim- 
ulations of galaxy mergers to show that there is indeed a 
better choice of dimensional scaling constants, at least for 
the set of simulations under consideration. Unfortunately 
our preferred choice of scaling constants are not directly 
observationally accessible, so in Section [5] we work out the 
relationship between our preferred definition of the scaling 
constants and the conventional choices for plausible bary- 
onic and dark matter distributions. Finally, in Section [6] we 
summarize our conclusions. 



2 THEORETICAL PRELIMINARIES 

When discussing single-component models, it is fairly clear 
how to make the definition of homology precise. Such a 
galaxy is fully described by its distribution function f(x, v) 
giving the density of stars with the given position x and 
velocity v. A set of galaxies is homologous if there exists 
a 'master' distribution function F where the distribution 
function of a specific galaxy / is obtainable from F via: 



M 

Wv 



3 F(x/R,v/V) 



(1) 



where R, V, and M are arbitrary scale constants, different 
for each galaxy. If one knows the universal distribution func- 
tion, three numbers suffice to completely describe a galaxy. 
This will be called scaling homology. 

Scaling homology is closely related to systems of units. 
It means that all galaxies are the same if one works with 
the correct system of units, where the units are allowed to 
change from galaxy to galaxy. 

It is very important to recognize at this point that these 
scaling constants are not necessarily equal to the half-mass 
radius, the velocity dispersion, and the total galaxy mass. 
For single component galaxies, these choices for the scaling 
constants are essentially the only reasonable ones: if galax- 
ies are homologous, then other choices will end up being 
simple functions of the half-mass radius, etc. However, for 
multi-component galaxies, our job is to take a set of galaxies 
and infer the definitions for R, V, and M that make all of 
their distribution functions derivable from a single 'master' 
distribution function. If no such definition exists, then the 
galaxies are not scaling-homologous. Section [4] is devoted 
to carrying out this analysis for a set of simulated galaxy 
merger remnants. 

Perhaps the most common method of defining homology 
is to define structure constants relating the true values of a 
galaxy's mass and luminosity to 'virial estimates' obtained 



L = cil e «r ef{ 



(2) 
(3) 



where L is the galaxy's luminosity, M is the mass, r e fl is the 
projected half-light radius, / e fi is the mean surface bright- 
ness inside r e g, ao is the projected aperture velocity disper- 
sion, and ci and C2 are structure constants defined by these 
two equations. The important point is that M is the true 
mass of the galaxy while Co and r e g are scale factors com- 
bined in such a way as to dimensionally produce a mass (up 
to physical constants). Therefore variation of C2 from galaxy 
to galaxy indicates that mass profiles of the two galaxies 
differ to produce a different relationship between the virial 
estimator and the actual mass. A set of galaxies is said to 
be homologous if these structure constants are the same for 
all galaxies. We denote this definition structure- constant ho- 
mology. 

It is easy to show that for single-component sys- 
tems, structure-constant homology and scaling homology 
are equivalent. Putting the two definitions into the same 
language, we have: 



P( x ) = / f(x,v)d v 

the total kinetic energy is 

KE=~ / \v\ 2 f(x,v)d 3 xd 3 v 

and the total potential energy is: 
pE = _G f p(x)p(y)d 3 xd 3 y 



(4) 



(5) 



- / \sb - y\ 

G f f{x,v)f{y,w)d 3 xd 3 yd 3 vd 3 w 

\ x - v\ 



(6) 



where G is Newton's gravitational constant. The scalar virial 
theorem is then: 



\v\ f(x, v) d xd v 

G f f(x,v)f(y,w)d 3 xd 3 yd 3 vd 3 w 
2 J \x-y\ 



(7) 



Let us define three arbitrary scaling constants R, V, and M, 
and dimensionless coordinates x — x/R and v = v/V. 



(8) 



J 

GR 5 V 6 f f(x, v)f(y, w) d 3 x d 3 y d 3 v d 3 w 
\x - y\ 



Using the definition of the dimensionless distribution func- 
tion in Equation [l] the virial theorem becomes: 



MV 2 I \v\ 2 F(x,v)d 3 xd 3 v 



GM 



F(x, v)F(y, w)d 3 x d 3 yd 3 v d 3 w 



using scaling constants. From Bender et al. ( 1992 \ 



2R J \x-y\ 
Therefore we define the structure constants 

\v\ 2 F(x , v) d 3 x d 3 v 



\v/V\ 2 f(x,v)d 3 xd 3 



(9) 



M 



(10) 
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and 



P = 



F(x, v)F{y, w) d x d 3 y d v d w 

I* - y\ 

R f f(x,v)f(y,w)d 3 xd 3 yd 3 vd 3 w 
M J \x-y\ 

With these definitions, the virial theorem is simply: 

aV = B — — 
M 2R 

and comparing to equation [3] reveals that: 
2a \ /Af bar \ / R \ (V^ 2 



f'2 



- I (' 

IB V 



(11) 



(12) 



(13) 



GB J \ M J \r cB J \a 
where Mbar refers to the total baryonic mass of the galaxy. 

2.1 Two Component Galaxies 

Describing a galaxy that contains both stellar and dark mat- 
ter requires specifying two distribution functions, f and g, 
with corresponding dimensionless functions F and G. 

Why do we use separate distribution functions for bary- 
onic matter and dark matter, but we do not separate the dif- 
ferent components of baryonic matter into young stars, old 
stars, disc stars, etc? The different stellar components have 
different formation histories and the moments of the distri- 
bution function that appear in the Jeans equations will be 
different for the different stellar components. Here we at- 
tempt to strike balance between accuracy and simplicity: 
our description of the system should be complex enough to 
be accurate, but simple enough to be informative. Hence we 
separate baryonic matter from dark matter because the two 
types of matter obey different physics and have very differ- 
ent histories. By including all baryonic matter in a single 
distribution function, we sacrifice some accuracy in order to 
make the description simple enough to be easily understood 
and manipulated. 

Allowing for a second component, the virial theorem 
takes the form: 



2RV 2 
GM 



\v\ 2 F(x, v) d 3 x d 3 v 



\v\ 2 G{x, v) d 3 x d 3 v 



F(x,v)F(y,w) 3 _ 3 _ 3 _ 3 . 
— - — ; — - — t-, d xd yd v d w 

\x-y\ 

, F(x,v)G(y,w) ,o o ,<>_ ,<» _ 
+2 / v \J 1, — '- d 3 xd 3 yd 3 vd A w 

\ x -y\ 

. G(x i)G(y w) , / W ^ /WV 
\x-y\ 
We define 

F(x,v)G{y,w)d 3 xd 3 yd 3 vd?w 

* - y\ 

2R f f(x,v)g(y,w)d 3 xd 3 yd 3 vd 3 w 
M J ' \x~~y~\ 

so that at last we have the virial theorem in the form: 



(14) 



(a f + a g )V 2 = (B f + B a+1 ) 



GM 
~2R 



(15) 



(16) 



From these expressions it is clear that for multi- 
component galaxies: 

• Scaling homology implies structure-constant homology. 



• Structure-constant homology does not imply scaling- 
homology because structure-constant homology only re- 
quires a small number of constraints on the definite integrals 
defining a, B, and 7. 

• A set of galaxies being structure-constant homol- 
ogous does not imply that each component is sepa- 
rately structure-constant homologous. Structure-constant- 
homologous galaxies can be built from non-homologous com- 
ponents if the change in a/ compensates for the change in 
a g and similarly for B. 

• The components of galaxies being structure-constant- 
homologous does not imply that the galaxies are structure- 
constant-homologous. 

• The structure constants for multi-component galaxies 
are not simply the sum of the structure constants of the 
individual components owing to the cross term 7. 

Finally, similar considerations for scaling-homology 
yield: 

• Scaling-homology of galaxies implies scaling-homology 
of the components. 

• Scaling-homology of the components does not imply 
scaling-homology of galaxies. If a set of galaxies is com- 
posed of scaling-homologous components, then the set of 
multi-component galaxies will be scaling-homologous only if 
the ratios of the scaling constants for each component are 
the same for all of the galaxies in the set. That is, Rf/R g 
must be the same for all of the galaxies. 

One may wonder whether the dark and luminous com- 
ponents should be allowed their own scaling constants. This 
idea is problematic because scaling the components sep- 
arately disturbs the dynamical equilibrium of the galaxy. 
Consider a galaxy in dynamical equilibrium where the bary- 
onic scale radius is allowed to shrink to zero, making the 
baryons into a point mass. The dark matter far outside the 
original baryonic scale radius will not be significantly af- 
fected by this operation. However, the dark matter near the 
centre (inside the original baryonic scale radius) will have 
more mass enclosed after the change. The new galaxy will 
obviously be out of dynamical equilibrium. 

We have seen from the above considerations that scaling 
homology is a stronger concept that structure-constant ho- 
mology in the sense that scaling homology implies structure- 
constant homology, but not vice versa. For two-component 
galaxies, the two components can undergo compensating 
changes to leave the structure constants unchanged. This 
does not seem to be a desirable property for a concept like 
homology. Finally, there is a certain theoretical elegance to 
the viewing homology as closely connected to the choice of 
the system of units — a set of galaxies is scaling homologous 
if all of the galaxies are the same when using the correct 
system of units. For all of these reasons we advocate that 
homology be taken to mean scaling homology. Structure- 
constant homology should be clearly denoted as a different 
concept. 



2.2 Tilt of the Fundamental Plane 

There are three independent ways to generate the tilt in the 
fundamental plane. Any or all of them may be at work. 

The first is a systematic change in the stellar mass- 
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to-light ratio with the mass of the host galaxy via trends 
in the age, metallicity, or initial mass function of the stellar 
population. In this case the stellar and dark matter distribu- 
tion functions are scaling (and therefore structure-constant) 
homologous because the masses of stars and dark matter 
components are unaffected. 

The second is a systematic change in the dark matter 
fraction near the centre of the galaxy. This is a mild form 
of non-homology which is achieved by changing the ratio 
of the mass scaling constants M^^/Mbm with galaxy mass, 
where Mdm is the total dark matter mass of the galaxy. 
Overall scaling- and structure-constant homology are vio- 
lated, but the baryonic components by themselves may be 
nearly scaling- and structure-constant-homologous. Further- 
more, as long as the stellar and dark matter remain nearly 
decoupled due to their very different half mass radii, then 
changing the dark matter fraction should not greatly dis- 
turb the dynamical equilibrium of the stellar material. In 
this case the stellar population mass-to-light ratio is identi- 
cal for all galaxies, but the total mass-to-light (both for the 
entire galaxy and within any radius) ratio changes system- 
atically with mass. 

The third possibility is that homology is flagrantly vi- 
olated and the baryonic components of galaxies are neither 
scaling nor structure-constant homologous. Even if the stel- 
lar population mass-to-light ratio is always the same, there 
is still total freedom in the total mass-to-light ratio owing 
to the total freedom in the choice of dark matter profile for 
each galaxy. It may change or be constant as a function of 
radius, and the total mass-to-light ratio may change or be 
constant for the entire set of galaxies. 

A great source of confusion is that the first and second 
possibilities are very different theoretically but rather simi- 
lar observationally. Everyone agrees that the first possibility 
involves homology and the third possibility involves non- 
homology. However, the second possibility is observationally 
very similar to the first possibility but theoretically very sim- 
ilar to the third possibility. This unfortunate circumstance 
makes it difficult to come to firm conclusions. 

2.3 Variable dark matter fraction and equilibrium 

There are two contributions to a total mass-to-light ratio: 
one from the characteristics of the stellar population and 
one due to the accompanying dark matter. If the mass-to- 
light ratio changes with galaxy mass because of a changing 
stellar population, then the FP will be tilted but the mass 
distribution may still be homologous. However, if the mass- 
to-light ratio changes due to changing dark matter fraction, 
then the galaxies must be non-homologous. 

Consider the spherically symmetric Jeans equation (e.g. 
Binney fc Tremaine|1987| >: 

d In v d In of 1 d<E> „ ,„„, 

-^ + -^ + 2^+^— =0 17 
dmr dmr cr; alnr 

where v is the number density of tracer particles, oy is the 
radial component of the velocity dispersion, $ is the total 
gravitational potential, /3 a = 1 — a% /of is the anisotropy pa- 
rameter, and oe is the tangential component of the velocity 
dispersion. 

The matter density is provided by two components, pbar 
and pdm for baryons and dark matter. The Poisson equation 



is linear so the solutions obey the superposition principle. 
Putting this into the Jeans equation results in: 

^ + ^ +2ft + J_(^I + « =0 (18) 
dmr dmr a* \ dmr dmr ) 

where $bar and <3>dm refer to the potential due to the 
baryonic and dark matter respectively and the constant b 
parametrizes the dark matter fraction. If a galaxy is in equi- 
librium for any value of b, then any change to b must be 
compensated by a change in the anisotropy /3 a or a change 
in the radial fall-off of density or velocity dispersion. 

Variations in the dark matter fraction and non- 
homology are the same concept. If the tilt in the FP is not 
due to systematic changes in the stellar population mass-to- 
light ratio, then it must be due to non-homology. There is 
no separate case where 'only' the dark matter fraction and 
therefore the observed mass-to-light ratios are changing. 

2.4 What constitutes non-homology? 

Given the somewhat confused use of the word homology in 
the present literature, it is perhaps easier to focus on what 
non-homology means. As discussed above, one way galaxies 
could be non-homologous is by having baryonic profiles that 
change shape from galaxy to galaxy. The fact that not all 
galaxies have the same Sersic index already tells us that 
the population of galaxies cannot be perfectly homologous. 
If the dark and luminous components scale separately from 
each other so that different galaxies have different amounts 
of dark matter within some radius, then the galaxies must 
also be considered non-homologous. 

Finally, recall that since structure-constant homology 
only places a few constraints on the definite integrals defin- 
ing a, P, and 7, the luminous and dark components of a set 
of galaxies can undergo compensating changes so that the 
galaxies would be deemed structure-constant homologous, 
but would not be scaling homologous. We use this observa- 
tion to argue that scaling homology is the more useful con- 
cept and the concept of structure-constant homology should 
be avoided. 



3 RECENT WORK 

Armed with a precise definition of the different flavors of 
homology we highlight a few recent observational and theo- 
retical results. We wish to give a brief overview of the state 
of the field, but unfortunately none of these studies will be 
directly relevant to the present one because they all use the 
half-light radius as the relevant length scale. In section|4]we 
will show that for a large set of hydrodynamic simulations, 
we obtain more instructive results by using r cq , the radius 
where the enclosed dark matter and luminous matter are 
equal. Simulated merger remnants look significantly more 
homologous when using r cq than when using r e ff. 

Thus far we have been concerned only with theoreti- 
cal galaxies: all quantities of interest can be measured to 
arbitrary precision. However, the only observationally ac- 
cessible quantities are projected ones. The measured mean 
surface brightness within the effective radius involves an in- 
tegral along the light of sight and therefore implicitly mixes 
information about the entire mass/light profile at all radii. 
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Furthermore, both dark matter and stellar populations con- 
tribute to the measured mass-to-light ratios of galaxies. 
Variations of mass-to-light ratios due to one or the other 
of these possibilities are very different theoretically, but ob- 
servationally similar. 

One may ask if we could take the results of an observa- 
tional or theoretical study discussed below, take a guess at 
the scaling of the density profile between r eq and r e ff, and 
then directly compare the results of previous studies to our 
own. Unfortunately, carrying out such a project would tell 
us something about our assumed scaling between r e fi and 
r cq , but it would not tell whether other authors would reach 
different conclusions if they use r cq instead of r e ff. In or- 
der to answer the latter question, we need access to the full 
density profiles as a function of radius in the case of numer- 
ical studies. We unfortunately must wait for other groups to 
analyse their simulations in terms of r eq in order to see if 
our findings hold true beyond our present set of simulations. 



3.1 Observational Studies 

There has been some success in recent years indicating that 
stellar population models have progressed to the point where 
it is possible to measure stellar mass-to-light ratios from 
spectra (e.g. |Graves fc Fa bcr 2010|. The usual assumption 
is that the initial mass function is universal, but recent ob- 



servations have cast doubt on that idea (van Dokkum & 
|Conroy|[2011| [Cappellari et al.|2012| . Furthermore, in spite 
of the many codes available to estimate stellar masses (e.g 



Bell et al.|2003||Kauffmann et al.|2 003 Gal lazzi et al.|2005| 



Blanton fc Roweis|2007^ " many are based on a rather small 
number of libraries of stellar spectra (e.g. |Bruzual fc Char~ 
lot|2003[ ). Thus one may be concerned about covariance be- 
tween the different estimates of stellar mass because their 
fundamental assumptions are similar. If this is the case, then 
agreement between the different estimates indicates that 
their assumptions are similar, not that the stellar masses 
have actually been accurately measured. 



Cappellari et al. ( 2006 1 used Schwarzschild ( 1979 \ mod- 



eling of SAURON ( |Bacon et al.|2001[ ) early-type galaxies to 
show that the so-called virial mass estimator (5r e ffO" 2 /G) is 
an excellent estimator of the dynamical mass within one ef- 
fective radius. The SAURON survey also has stellar spectra 
over the whole galaxy image out to approximately one ef- 
fective radius, and is thus excellently placed to use stellar 
population synthesis to finally break the degeneracy between 
dark matter fraction and stellar population mass-to-light ra- 
tio. Their Figure 17 unfortunately does not permit a single 
firm conclusion. For dynamical mass-to-light ratios below 
2.8, there is a linear relationship between the stellar mass- 
to-light ratio and the dynamical mass-to-light ratio, indi- 
cating that the dark matter fraction does not change. For 
larger mass-to-light ratios, there is essentially a single stellar 
mass-to-light ratio, independent of the dynamical mass-to- 
light ratio, indicating that either the initial mass function 
or the the dark matter fraction does change. 



Bolton et al. ( 2008 1 and Koopmans et al. ( 2009 ) used 



strong gravitational lenses to argue that if one replaces the 
luminosity surface density in the FP with the mass surface 
density, one gets an untilted fundamental mass plane that 
is in full agreement with the virial expectation. Therefore 
the tilt in the FP must be due to variations in mass-to-light 



ratio. However, there is little guidance from strong gravi- 
tational lensing about whether the total mass-to-light vari- 
ation is due to changing dark matter fractions, changing 
stellar populations, or a changing initial mass function. 



Taylor et al. (20101 used SDSS data to compare stel- 



lar mass and dynamical mass estimates for a set of galax- 
ies. They found that the two were related, but non-linearly, 
and that the non-linearity could be removed by defining a 
'structure-corrected' dynamical mass that takes the Sersic 
parameter n of the galactic light into account. Thus they 
conclude that elliptical galaxies are not homologous. Taylor 
et al. (20101 made the conventional choice of scaling dif- 



ferent galaxies by their half-light radii and central velocity 
dispersion in order to evaluate the question of whether or 
not they are homologous. 



Graves & Faber (20101 study stacked SDSS spectra of 



early- type galaxies as a function of their position in the Fun- 
damental Plane. They found that the variation in the dy- 
namical mass-to-light ratios was too large to be explained 
by variations in the stellar mass-to-light ratios alone. This 
indicates varying dark matter fractions and is interpreted as 
non-homology in the galaxy population. 



Auger et al. ( 2010 1 used multi-band Hubble Space Tele- 



scope images to construct stellar masses from stellar popu- 
lation models for 73 massive early-type galaxies identified 
as strong gravitational lenses as part of the SLACS survey 



(Czoske et al. 2008). Combining this information with ve- 



locity dispersions from the Sloan Digital Sky Survey (SDSS) 
(York et al. 20001, they concluded that the central dark 



matter fractions of massive early-type galaxies increase sys- 
tematically with galaxy mass if one assumes a universal ini- 
tial mass function. They found that the stellar mass-to-light 
ratio was essentially constant with galaxy mass. |Barnabe| 
et al. (20111 carried out dynamical modeling of a subset 



of these galaxies using ground-based integral-field-unit data 
and reached the same conclusion. 



Grillo (20101 analyzed SDSS velocity dispersions along 



with stellar masses obtained from the JHU/MPA value- 
added galaxy catalog]]] Dynamical considerations led them 
to conclude that the ratio of dynamical mass to dark mass 
was constant and they argued that the tilt of the FP is 
entirely due to stellar population effects. |Grillo fc Gobat| 
( |2010[ ) used Schwarzschild ( 1979 \ models of 13 bright galax- 



ies in the Coma cluster constructed by Thomas et al. ( 2007 
|2009[ ) along with stellar masses from SDSS to reach the same 
conclusion. 

After many years of effort, it seems that stellar pop- 
ulation modeling and dynamical modeling have progressed 
to the point where it is possible to separate the effects of 
changing dark matter fractions from that of changing stellar 
populations. However, recent studies have concluded both 
that galaxies are homologous and that they are not homol- 
ogous. In addition, most applications of stellar population 
modeling require the assumption of a universal initial mass 



function. However, work by van Dokkum & Conroy (20111 



and Cappellari et al. (20121 have cast doubt on this long- 



held assumption. If the initial mass function is truly not 
universal, both dynamical modeling and stellar population 
modeling will be more difficult. 



http: / /www. mpa-garching.mpg.de/SDSS/ 
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3.2 Theoretical Studies 



Gonzalez-Garcfa & van Albada (20031 and Boylan-Kolchin 



et al. (20051 studied dissipationless merger simulations to 



determine whether mergers leave the FP intact. Gonzalez- 
Garcfa & van Albada (2003) started with structure-constant 



homologous systems and found that their remnants did not 
share the property. The merger remnants do lie on the FP, 
evidently because the two structure constants are able to un- 
dergo compensating changes since they only appear in a mul- 
tiplicative combination in the FP relation (their Equation 
6). Boylan-Kolchin et al. (20051 studied a smaller number 
of much higher resolution simulations and found that both 
structure-constant homology and the FP were preserved by 
the merger process at the 10% level. They did, however, find 
a significant violation of homology for radial-orbit mergers. 



Dekel & Cox ( 2006 1 used numerical simulations to study 



dissipation as the potential origin of the FP tilt. They as- 
sumed that all observable galaxy properties are power-law 
functions of the mass of the galaxy and the gas fraction, 
and then used several merger simulations with different gas 
fractions and masses in order to fix the power-law exponents 
and constants in their model. They were left with a set of 
constraints on the power-law exponents such that for any 
observed FP tilt, one could solve for the properties of the 
progenitor galaxies (e.g., gas fraction) as a function of mass. 
This of course depends on their merger simulations being 
good representatives of all mergers of that mass and gas 
fraction. 

|Robertson et alT] ( |2006| ) used a large number of dissi- 
pational simulations of galaxy mergers, assumed a constant 
stellar mass-to-light ratio and concluded that with a suffi- 
cient amount of gaseous dissipation, merger remnants nat- 
urally produce a tilted FP. Since they assumed a constant 
stellar mass-to-light ratio, this tilt must have been caused 
by structural non-homology/ varying dark matter fraction. 



Hopkins et al. (20081 performed an extensive study of 



dissipation and the origin of the fundamental plane. They 
seem to be concerned with structure-constant homology, but 
dismiss non-homology as the source of the FP tilt. They as- 
sert that gas dissipation causes a varying dark-matter frac- 
tion inside the effective radius, which in turn tilts the FP. As 
stated above, we note that a varying dark-matter fraction is 
not distinct from non-homology. However, it can perhaps be 
considered to be a weak form of non-homology where only 
one of many possible degrees of freedom is exercised. 



3.3 Relation to the present study 

It is conceivable that any or all of these studies would reach 
a different conclusion if they were to use our definition of 
length, velocity, and mass scales rather than the usual ob- 
servationally motivated definitions. In order to make this 
determination, we need more information than is available 
in the published papers, so we must wait for other groups 
to try analyzing their simulations and observations in terms 
of r cq rather than r of f . 

We hasten to add that we are not claiming that any 
previous authors have made any mistake. Other researchers 
have correctly formulated and answered questions about 
the scaling relations of galaxies. Given the assumption that 
the relevant length scale was the baryonic half-mass radius, 



these studies generally found measurable non-homology in 
the galaxy population that manifested itself by a changing 
central dark matter fraction. 

We are pointing out that the choice of r c ff as the rele- 
vant length scale is arbitrary, and other choices may be more 
instructive in understanding galaxies and galaxy merger 
remnants. On the bases of our simulations, we have found 
that scaling the simulations by r eq gives very instructive 
results. We hope that other researchers will find this termi- 
nology similarly useful. 



4 APPLICATION TO SIMULATIONS 

As already noted, the dimensional scaling constants used 
in the definition of scaling-homology are not required to be 
the conventional choices of half-mass radius, central veloc- 
ity dispersion, etc. In this section we consider a set of simu- 
lated galaxy merger remnants and search for a definition of 
the dimensional scaling constants so that the dimensionless 
properties of the set of remnants are as similar as possible. 

Once one has arrived at a definition for the scaling con- 
stants that seems to work well, one must compute the value 
for a given galaxy. Therefore we are restricting ourselves to 
simple, physically-motivated definitions of the dimensional 
scaling constants, such as the baryonic half-mass radius, the 
total half-mass radius, the radius where the enclosed bary- 
onic and dark masses are equal, etc. 

One could imagine carrying out a high-dimensional min- 
imization where the function to be minimized is the sum of 
the total deviations of the scaled density and velocity profiles 
from the mean scaled profiles for all galaxies. The free pa- 
rameters in the minimization are the three scaling constants 
for each galaxy. This would effectively tell you whether or 
not there exists a definition of the scaling constants that 
makes the galaxies scaling-homologous. However, it would 
not give a simple way to compute the dimensional scaling 
constants for a given galaxy. One would know the value of 
the scaling constants for a galaxy in the 'training' set: the 
fitting procedure would provide it. But one would not be 
able to compute the dimensional constant for a new galaxy, 
not in the training set. 

Using numerical simulations allows us to quickly explore 
many possible definitions of the scaling constants with little 
uncertainty in, for example, the total galaxy mass or velocity 
dispersion. The simulated galaxies are of course imperfect 
stand-ins for actual galaxies, but our hope is that working 
with simulated galaxies will allow us to eliminate many of 
the definitions that are possible and focus on one (or a few) 
that seem to work well. Only the most successful choices 
when compared to simulations will be considered including 
the additional complications inherent to observational data. 

There are two sets of simulations of binary mergers con- 
sidered here. The Sbc series consists of equal mass mergers 
between two identical gas-rich progenitors on a variety of or- 
bits. The G series consists of mergers of galaxies with lower 
gas fractions and a wide range mass ratios using only a few 



different orbits. Cox et al. (20061 contains a full description 



of the simulations. 

We considered eight possible radius scales including 
stellar, dark, and total 3D effective radii and projected ef- 
fective radii, the radius where stellar mass density and dark 
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mass density are equal, and finally the radius where the en- 
closed stellar mass and enclosed dark mass are equal. We 
evaluated thirty mass scales including the stellar, dark, and 
total mass, the total of each type of mass within each of 
four radii (stellar effective radius, dark effective radius, ra- 
dius where stellar and dark densities are equal, and the ra- 
dius where the enclosed stellar and dark masses are equal), 
and the mass scale that sets each of the three types of mass 
density to one at each of the four radii. The last dimensional 
choice is either a velocity or a time, and we tried the cen- 
tral 3D velocity dispersion, the projected velocity dispersion, 
and dynamical times at each of the above listed radii. 

Evaluating one 'model' requires a choice of one of each 
of the three dimensional constants: length, mass, and veloc- 
ity/time. This leads to a very large number of combinations. 
Fortunately, only a few choices work better than the conven- 
tional observationally motivated choices of baryonic half- 
mass radius, central velocity dispersion, and total baryonic 
mass. In this paper we plot only the conventional choices 
and the best alternate definition we found: the length scale 
is r oq , the radius such that the enclosed dark and baryonic 



^dyn — 



so that 



3tt 
32Gp 



8GM C , 



(19) 



(20) 



To evaluate each choice of the scale constants, we plot 
dimensionless density and velocity dispersion profiles for all 
galaxies: 

R 3 p(r/R) 



p(r) 



M 

°(r/R) 
V 



(21) 



(22) 



where r is the magnitude of the radius vector. 

For simplicity we use spherical apertures to convert the 
3D density and velocity profiles of the simulated galaxies 
into the ID profiles shown here. Using ellipsoidal profiles 
instead does not change our conclusions. The velocity dis- 
persion and rotation information is given in spherical coor- 
dinates r, 0, and <f> where <j> is the azimuthal angle and the 
coordinate system is defined by the shape of the stellar rem- 
nant so that the z axis corresponds to the short axis of the 
remnant. 

If the simulated galaxies are perfectly scaling- 
homologous and we have hit upon the correct definition of 
the dimensional scaling constants, then the dimensionless 
density and velocity dispersion profiles will be exactly the 
same for all galaxies. Non-scaling-homologous galaxies will 
lead to large variations in the profiles from galaxy to galaxy. 
Therefore we plot the dimensionless profiles for all galaxies 
(with some transparency in each line so that the effects of 
over-plotting are minimized) and seek the definition of scal- 
ing constants that leads to the least variation from galaxy 
to galaxy. 

For each simulation, we must estimate various observa- 
tionally relevant quantities from simulation snapshots. For 



real galaxies, observables such as the half-light radius, aper- 
ture velocity dispersion, and luminosity are affected by star 
formation history, age gradients, metallicity gradients, and 
dust. We ignore the majority of these complications and 
make simple assumptions that allow us to focus on the es- 
sential physics that determine the structure of merger rem- 
nants. In the future, a more detailed study using a code that 
allows us to take the effects stellar populations and dust into 



account (e.g. Jonsson 2006 Jonsson et al. 2010 1 would be 
worthwhile. 

For the moment, we assume a constant stellar mass-to- 
light ratio and no dust scattering or absorption allowing for 
a simple relationship between mass and light. To estimate 
3D and aperture velocity dispersions, we perform a maxi- 
mum likelihood fit of a Gaussian distribution to the simula- 
tion particle velocities. The stellar remnants are very close 
to oblate spheroids, so to estimate the projected half-light 
radius we use a "typical" viewing angle 45 degrees between 
the pole-on and face-on views. We fit an ellipse to the pro- 
jected iso-density contours and use the geometric mean of 
the semi-major and semi-minor axes as the half-light radius. 

Figure[T]shows the dimensionless density profiles for the 
usual choice of the stellar effective radius for the length scale 
and the total stellar mass for the mass scale. 

Figure [2] shows a different choice of length and mass 
scales. The length scale is r eq , the radius at which the en- 
closed stellar and dark masses are equal, and the mass scale 
is the stellar mass enclosed within r cq . In comparing to 
the mass profiles shown in Figure [I] we must ensure that 
the same physical radial range of the galaxies are probed. 
Multiplying the chosen units by a factor common to all 
galaxies does not change the result, so we normalize to the 
5'6c201a-u4 simulation. For this simulation, r e g = 0.437r cq 



and M ba 



1.54M e , 



The simulated galaxy merger remnants analysed here 
are not perfectly scaling-homologous. This is clear as soon 
as one realizes that the central slopes of the dark matter 
profiles are different between the G and Sbc series of simu- 
lations. Changing the scale radius can only shift curves from 
side-to-side in Figures [I] and [2] it cannot change the slopes. 
Even so, the radius and mass scale choice made in Figure [2] 
represents an impressively high degree of regularity among 
galaxy merger remnants given the variation in mass ratios, 
orbits, and gas fractions. 

The total mass density profile for all of the simulations 
is very close to an isothermal p oc r~ 2 profile. This is in good 
agreement with recent studies of the mass profiles of ellip- 
tical galaxies using strong and weak gravitational lowing 
( |Gavazzi et al.||2007| |2008 | |Bolton et aT1|2008 |l, dynamical 



studies of the line-of-sight velocity profiles of giant elliptical 
galaxies (Gerhard et al. 20011, and studies of the density 



and temperature profiles of x-ray gas in elliptical galaxies 
spanning a wide range in mass ( [Humphrey fc Buote|2010| . 

Table [I] lists the RMS variation among all simulations 
for each physical quantity considered in this section. We take 
the value of the non-dimensional quantity under considera- 
tion at the appropriate scale radius (either r e ff or 0.437r cq ) 
and compute the RMS variation of the resulting set of num- 
bers. This serves as a indication of the galaxy-to-galaxy vari- 
ation of, for example, the density profiles after applying a 
given set of scaling constants, and therefore of the utility of 
that particular set of scaling constants. 
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Figure 1. The dimensionlcss mass density profiles for all G and 
Sbc series simulations using the typical observationally motivated 
choices of stellar half mass radius and stellar mass for the length 
and mass scales. Blue lines are for Sbc simulations (with gas rich 
progenitors) while red lines are G series simulations. From the 
top, the total mass density, the stellar mass density, and the dark 
matter mass density. All mass profiles have been multiplied by r 2 
to flatten them out. The stellar components of all of the simula- 
tions are strikingly similar. However, the chosen mass and length 
scales have nothing to do with the dark matter, so there is sig- 
nificant variation among the dark matter profiles of the different 
simulations. There is also a systematic change in the slope of the 
dark matter profiles between Sbc and G series simulations. No 
change of units will remove this difference. 



Figure 2. The dimensionlcss mass density profiles for all G and 
Sbc series simulations using our preferred choice of radius scale 
of 0.437r cq , where r eq is the 3D radius where the enclosed stellar 
mass and dark mass are equal, and the mass scale is the 1.54 
times the stellar mass enclosed within r eq . The factors of 0.437 
and 1.54 are chosen to ensure that these plots cover nearly the 
same physical region of the galaxy as Figure [T] With this choice 
of length and mass scale, the merger remnants show remarkable 
regularity. The dark matter mass profiles show 6% root-mean- 
square variation near the half-mass radius and the baryonic mass 
profiles show 17.7% RMS variation (compared to 51% and 18.5% 
respectively for the conventional observationally motivated set of 
scaling constants) in spite of the wide range of masses, gas frac- 
tions, orbits, and feedback recipes employed in these simulations. 



Figure [3] shows remnant stellar velocity profiles for the 
typical choice of scaling constants: the stellar half-mass ra- 
dius and the central projected velocity dispersion within 
r e ff/8. Figure [4] shows stellar velocity profiles when scaled 
by r cq and the dynamical time at r eq . 

The velocity dispersion anisotropy profiles provide a 
convenient limit on the extent to which the velocity struc- 
ture of the simulated remnants can be brought into a com- 
mon system. The anisotropy is dimensionless, hence the 



only choice to be made is the length scaling constant. The 
anisotropy /3 a is usually defined for a spherically symmetric 
object. In that case, ag = cr^ and v r — vg = — The 
merger remnants presented here are not spherically symmet- 
ric so we define 

= 1 - (23) 
and 

04 = 1- 4/ ff r (24) 
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Figure 3. Dimensionless stellar velocity profiles with scale ra- 
dius and velocity taken to be the usual observationally motivated 
ones: ro is the is the stellar half-mass radius and cto is the pro- 
jected velocity dispersion within an aperture of r e ff/8. From the 
top: 3D radial velocity dispersion, 3D velocity dispersion in the 
direction, 3D velocity dispersion in the azimuthal direction, and 
mean streaming velocity in the azimuthal direction. There is sig- 
nificant diversity in the velocity structure of simulated galaxy 
merger remnants. 



Figure 4. Dimensionless stellar velocity profiles with scale radius 
and time taken to be our preferred choices of r oq and t cq . This is 
an improvement over the conventional choices for scaling radius 
and velocity shown in Figure [3] the root-mean-square variation 
of ay, for example, at the scale radius goes from 27% to 15% in 
adopting the new set of scaling constants. The RMS variations of 
all quantities shown here are listed in Table [l] 
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Table 1. Root-mean-square variation of non-dimensional quan- 
tities among all simulations at the given scale radius. The con- 
ventional observationally-motivatcd choice of scaling constants is 
used for the column on the left (r e g, Mt, ar , and (Tq), while our 
preferred choice of scaling constants based on the radius at which 
the enclosed dark and luminous masses are equal (r eq , M eq , and 
t eq ) is used for the column on the right. The quantities listed are 
the total mass density, the baryonic mass density, the dark matter 
mass density, the radial 3D stellar velocity dispersion, the 3D stel- 
lar velocity dispersion in the azimuthal direction, the 3D stellar 
velocity dispersion in the 8 direction, the mean streaming veloc- 
ity in the azimuthal direction, the specific kinetic energy of the 
stellar component, the specific kinetic energy of the dark matter 
component, the stellar velocity anisotropy in the azimuthal direc- 
tion, and the stellar velocity anisotropy in the 8 direction. All of 
the values are expressed as a percentage of the mean value for all 
galaxies except for the anisotropy /3, where we do not divide by 
the mean because anisotropy is not positive definite. 



Quantity 


r e ff based scaling 


r eq based scaling 


Pbar+PDM 


30.1% 


10.1% 


Pbar 


18.5% 


17.7% 


PDM 


51.3% 


6.04% 


° V 


26.9% 


15.0% 


0~tj, 


37.6% 


21.0% 


&e 


40.8% 


19.8% 




188% 


176% 


Ek/m (stars) 


94.9% 


11.0% 


Ek/m (dark matter) 


64.6% 


14.1% 


A* 


0.263 


0.238 


A? 


0.176 


0.166 



Figure [5] shows both /3$ and j3$. 

Virtually all of the remnants show stellar anisotropy 
profiles that rise with radius, in some cases to values near 1 
(completely radial velocity dispersion). This is easily under- 
stood in terms of the dynamical origin of the stars at large 
radius in merger remnants ( jDekel et al. 2005 1 . Stars that end 
up at large radius are stars that were flung out from near the 
centre of mass when the merging galaxies underwent their 
first close pass. They have low angular momentum because 
they started out within approximately one effective radius, 
hence the anisotropy measured at large radius will be very 
radial. 

The slopes of the anisotropy profiles are similar from 
galaxy to galaxy, but they are shallow. Changing the scale 
radius slides the curves horizontally, so any scale radius that 
would bring all the curves into agreement would be forced 
to vary dramatically from galaxy to galaxy. 

We have seen that the density profiles of simulated rem- 
nants are remarkably consistent, but the velocity profiles 
show significant variation. Figure [5] suggests that the rea- 
son for the variation in the velocity profiles is the direction 
of the velocities of the stars. That is, there is a common 
kinetic energy profile, but the velocity profiles are different 
because some remnants are radially anisotropic, some are 
tangentially anisotropic, and some are rotating. 

The specific kinetic energy is: 




Ek/m = -(«0 + a r + a B + <r ) 



(25) 



where we have assumed Gaussian velocity distributions and 



Figure 5. Stellar velocity dispersion anisotropy profiles for the 
usual choice of scale radius: the baryonic half-mass radius. This 
shows the reason why there is diversity in the velocity profiles of 
galaxy remnants in spite of the regularity in the density profiles. 
The virial theorem constrains the total kinetic energy of the sys- 
tem, but does not determine the direction of the velocities of the 
stars. Simulated remnants make use of this freedom and have a 
wide range of anisotropy profiles. Since the velocity anisotropy 
is dimensionless, changing dimensional scaling constants for each 
galaxy only moves each line horizontally. The shallow slope of 
the lines means that changing the definition of dimensional scal- 
ing constants cannot significantly change this plot to make the 
merger remnants more homologous. 



neglected v r and ve because the remnants are in steady- 
state. 

Figure [6] shows kinetic energy profiles using the stellar 
half-mass radius and central projected velocity dispersion 
as scale constants. A common kinetic energy profile is not 
apparent with this choice of scale constants. 

Figure [7] shows the dimensionless kinetic energy scaled 
by the mass-equality radius and the dynamical time at the 
mass-equality radius. The common energy structure of the 
merger remnants is apparent. In the case of the Sbc simu- 
lations there is as little as 3% RMS variation at the scale 
radius. There is again an offset between the two sets of sim- 
ulations but nevertheless the lack of variation between sim- 
ulations of a given set is remarkable. 

It is clear that dark matter plays a significant but sub- 
dominant role in the dynamical equilibrium of the luminous 
parts of galaxies, and that the baryonic and dark matter 
components of galaxies scale very differently with radius. 
Observational and theoretical studies to date have almost 
exclusively used the classic observationally accessible half- 
light radius as the relevant scaling parameter to compare 
galaxies of different sizes and masses. However, the value 
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Figure 6. Specific kinetic energy profiles with scale constants ro 
equal to the stellar half mass radius and to = to/hq. The upper 
panel shows the stellar component while the lower panel shows 
the dark matter component. There appears to be great diversity 
in the kinetic energy profiles of galaxy remnants. 

of this parameter involves only the luminous matter. It is 
therefore quite plausible that adopting a different scaling 
radius will lead to different conclusions about homology: a 
set of galaxies (real or simulated) may be closer to homology 
when using r cq as the scaling parameter. 

Unfortunately it is very difficult to determine what con- 
clusions other studies would have reached if they had used 
r eq instead of r c s without access to the full density and light 
profiles of the objects that were used in the study. One may 
be able to make reasonable extrapolations of the profiles 
from r c fr to r cq , but this would not test the hypothesis. It is 
the conclusion of the present study that sets of galaxies may 
look more homologous or less homologous depending on the 
scaling radius that one uses, so making an assumption about 
the mass profiles in order to translate results from one radius 
to another would not test the hypothesis. One needs access 
to the full mass profile (in the case of simulated galaxies) or 
light profile (in the case of observational studies). On this 
point we hope that other observers and theorists use the rel- 
atively simple procedure of scaling by r cq to see, first, if it 
makes their results easier to understand, and, second, if the 
results of our present study hold for broader sets of simula- 
tions and observations. We believe that adopting r cq rather 
than r e s would result in a significant gain in understanding. 



5 OBSERVATIONAL IMPLICATIONS 

As we have already mentioned, it is unfortunately the case 
that what we find to be the best definitions for the dimen- 
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Figure 7. Specific kinetic energy profiles with scale constants 
taken to be our preferred choices of r eq and the dynamical time at 
r eq . The upper panel shows the stellar component and the lower 
panel shows the dark matter component. These kinetic energy 
profiles show very consistent behaviour given the great diversity 
of the origins of the remnants. 

sional scaling constants that make simulated merger rem- 
nants close to homologous are not directly observable. In 
this section we define a plausible model for the density pro- 
files of baryonic and dark matter in a two component galaxy 
and find the relationship between our preferred dimensional 
scaling constants and the observationally accessible values. 
That is, we find r cq /r a g and Af cq /Afb a r for a range of plau- 
sible mass models. 

It is possible to plot values of r oq and M eq for our sim- 
ulated galaxy merger remnants but there are several prob- 
lems with this idea. Most importantly, anyone who wished 
to use these values of r eq and M eq to interpret observations 
would automatically be subject to all of the uncertainties 
and deficiencies of numerical simulations, for example in the 
star formation law to name one issue among many. Second, 
because our galaxies are so close to homologous when r cq 
is used as a scaling parameter, the values of r cq /r c ff and 
M oq /Mbar are nearly constant. Third, this particular set of 
simulated merger remnants do not span the full space of 
observationally relevant parameters, in the Sersic index for 
example. 

Therefore we seek a simple, easily understood, widely 
used, and observationally relevant mass model for which we 
can compute r eq and M cq . It is easy to understand the as- 
sumptions that enter such a model and, importantly, the 
model is not dependent on the details of our simulations 
even though the motivation to compute r eq and M eq arose 
from analysis of our simulations. 

It is a relatively simple matter to calculate r eq and M cq 
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for any given combination of dark matter and baryonic mass 
profiles. We provide plots and fitting formulas for the most 
common cases in the hope that these formulas may be used 
in both observational and theoretical efforts to understand 
the scaling relations of galaxies. It is the conclusion of this 
study that different choices of the scaling parameters can 
lead to different conclusions about homology, so this is po- 
tentially an important point for both theoretical and obser- 
vational studies of galaxy scaling relations. 

For the baryonic component of our mass model, we use 
deprojected Sersic profiles ( Sersic|1968 1: 



I(r p ) = Io exp 



m 



l/n 



(26) 



where I is the surface brightness, r v is the projected radius, 
7o is the central surface brightness, and r s is a scale radius 
chosen to ensure that the half-light radius is equal to some 
desired value. The Sersic profile is defined for the projected 
surface brightness of a galaxy image. Approximations to the 
three-dimensional luminosity density that lead to the Sersic 
profile under the assumption of spherical symmetry have 
been available for some time (Prugniel fc Simien|1997[ Lima 



Neto et al.|1999||Trujillo et al.|2002||Mamon fc Lokas|2005 l 

Exact analytic deprojections of the Sersic profile for integer 
values of the Sersic index n have also been available for some 



time (Mazure & Capelato 


20021. Recently Baes & Gentile 


( 201 ) have generalized the 


Mazure & Capelato|(|2002|) result 



for rational values of the Sersic index, a very significant step 
forward for practical applications of the formulas since ratio- 
nal numbers can approximate any desired value arbitrarily 
well. These involve special functions, but numerical imple- 
mentations of the necessary functions are available. |Baes fc| 
|Gentile| ( |2010[ > also provide formulas for the case where the 
Sersic index is real- valued rather than rational, but numeri- 
cal implementations of the required special functions are not 
readily available. We use the |Baes fc Gentile] ( |2010"[ ) formulas 
for exact deprojections of the Sersic profile. 

We allow two possibilities for the dark matter halo. The 



first is the familiar NFW formula ( Navarro et al.|1996 l: 

Pnfw 



p(r) 



x(l + x) 2 



(27) 



where p is the mass density, x = r/fNFW, ?"npw is the scale 
radius where the profile switches from a logarithmic slope of 
-I to -3, and Pnfw is a constant that sets the dark matter 
fraction of the galaxy as a function of radius. 

There is good evidence from strong and weak gravita- 
tional lensing studies that the total mass density of elliptical 
galaxies is remarkably close to that of a singular isothermal 
sphere (SIS) over a large range in radii ( |Gavazzi et al.|2007 



2008 Bolton et al.||2008). Therefore the second possibility 



for the dark matter halo is the mass density such that the 
total mass density is given by: 



f>r(r) = 



27rGr 2 



(28) 



where s parametrizes the depth of the potential well. 

Figure[8]gives r cq /r ff and M eq /Mbar for a range of mod- 
els assuming an NFW halo parametrized by the shape of the 
baryonic mass profile and the dark matter fraction within 
one baryonic effective radius. For a fixed baryonic mass and 
effective radius, the dark matter fraction within one bary- 
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Figure 8. The relationship between our preferred dimensional 
scaling constants and directly observable quantities for a vari- 
ety of mass models. The upper panel shows r eq /r e jr while the 
lower panel shows M oq /M E ff, both as a function of the 3D dark 
matter fraction within the projected baryonic half-light radius. 
The parameter n is defined in equation |26| and gives the shape 
of the baryonic density profile. The dark matter fraction within 
one baryonic effective radius changes as a result of changes to 
the NFW density parameter Pnfw f° r a fixed baryonic model. 
The dominant consideration is the dark matter fraction, while 
the shape of the baryonic mass profile has a smaller effect. 



onic projected half-light radius is adjusted by changing the 
NFW density parameter pnfw- The results depend on the 
ratio of the tnfw /r e s which indicates whether the density 
profile of the dark matter halo changes logarithmic slope 
near the baryonic component or far outside of it. The de- 
pendence is weak as long as rNFw/feff is sufficiently large, 
with the values of r cq and M eq approaching constant val- 
ues. The plots and fitting formulas presented here assume 
rNFwAcft = 20. The values of r eq are ~ 10% smaller for 
rNFwAcft = oo and ~ 10% larger for rNFw/?"cff = 10, 
with little dependence on Sersic index or central dark mat- 
ter fraction. When rNFw/fcff approaches unity, r eq depends 
sensitively on the ratio, but that corresponds to halo con- 
centrations C — r v i r i a i/rNFW approaching 100, far outside 



the range predicted by simulations for main haloes ( Bullock 
|et al .||2001[), although tidal stripping can lead to very large 



values of the concentration for sub-haloes (Diemand et al. 
|2008| ). 

The value of r eq for the case of an NFW halo may be 
approximated as: 



r e ff 



logic, .A(/dm) + .M/dm) lo gio( n / 3 ) 
+/ 3 (/D M )log? (n/3) (29) 
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Table 2. Fitting formula constants for NFW haloes. 



Name 


Value 


ki 


1.183 


A-2 


0.7252 


A3 


0.7241 


ki 


0.5476 


fcs 


2.752 


ke 


-0.1582 


k 7 


0.09714 


kg 


-0.6602 


kg 


0.7261 


kio 


-0.05721 


fen 


-0.005238 


ki2 


1.370 


ki3 


-1.319 



Table 3. Fitting formula constants for SIS haloes. 



Name 


Value for /dm 


< 0.55 


Value for /dm >= 0.55 


fcl 


2.546 




2.546 


A-2 


0.6363 




0.6363 


fcs 


0.1754 




0.1754 


A'4 


1.097 




1.097 


At, 


3.779 




3.779 


A(i 


-0.3725 




-43.72 


A 7 


0.2549 




38.30 


Ah 


-1.007 







An 


0.7579 







fcio 


-0.1034 




47.96 


An 


0.8789 




-42.32 


A'12 


-2.0134 







fcl3 


2.094 







h{x) 
h{x) 
h(x) 


= k 1 (x/k 2 ) k3 W- ki( - x/k2)k5 

= k 6 + k 7 (x/0.5) + fe 8 (x/0.5) 2 + k 9 (x/0.5) s 

= k w + feu(a:/0.5) + k 12 (x/0.5) 2 + k 13 (x/0.5) 3 



(30) 
(31) 



where n is the Sersic index, /dm is the 3D fraction of dark 
matter to total matter within the projected baryonic half- 
light radius, and the constants fci through fci3 are given in 
Table [2] This approximation is accurate to 8% for Sersic 
indices between 1 and 8 and and dark matter fractions 
between 5% and 70%. This formula gives the values for 
rNFw/fcft = 20 and is subject to modification described 
above. 

Figure [9] gives Teq/^efF and iVf eq / 'M\> aT for different cen- 
tral dark matter fractions under the assumption that the 
total density profile is that of a singular isothermal sphere. 
For a fixed baryonic mass and effective radius, the dark mat- 
ter fraction within one baryonic projected half-light radius 
is adjusted by changing the SIS parameter s. 

In this case r Gq and M Gq may be approximated by the 
same set of equations as for the case of an NFW halo: equa- 
tions [29] through [32] Constants for the equations are given 
in Tabic [3] This approximation is accurate to 5% for Sersic 
indices between 1 and 8 and dark matter fractions between 
0.05 and 0.7, except near the transition between the two 
sets of fitting formula constants at /dm = 0.55, where the 
approximation is only good to 20%. 



2 




1 




0.5 






— n=l 




n=2 


0.2 


— n=4 




n=6 





Figure 9. The relationship between our preferred dimensional 
scaling constants and directly observable quantities for a vari- 
ety of mass models. The upper panel shows r eq /r e g while the 
lower panel shows M cq /Af e ff , both as a function of the 3D dark 
matter fraction within the projected baryonic half-light radius. 
The parameter n is defined in equation [26] and gives the shape of 
the baryonic density profile. The dark matter fraction within one 
baryonic effective radius changes as a result of changes to the SIS 
density parameter s for a fixed baryonic model. The dominant 
consideration is the dark matter fraction, while the shape of the 
baryonic mass profile has a smaller effect. 



6 CONCLUSIONS 

We have considered the exact meaning of the concept of ho- 
mology in the context of multi-component galaxies. A par- 
ticularly useful formulation of the concept is termed scaling- 
homology. Using this definition, homology means that all 
galaxies are the same when scaled using the correct choice 
of length, mass, and time units, where the units are allowed 
to change from galaxy to galaxy. The three dimensional con- 
stants serve to fully specify a galaxy once the 'master' dis- 
tribution function that describes all galaxies is known. 

We used a set of hydrodynamical binary galaxy merger 
simulations spanning a range of masses, mass ratios, orbits, 
gas fractions, and supernova feedback prescriptions (but no 
AGN feedback) to determine the definition of dimensional 
scaling constants that served to make the remnants as close 
as possible to homologous. We found that the best defini- 
tion of length, mass, and time units are r oq , the radius within 
which dark and baryonic masses are equal, the mass within 
r cq , and the dynamical time measured at r cq . This is in 
contrast to the usual observationally motivated definitions 
of dimensional scaling constants: the baryonic half-light ra- 
dius, the baryonic mass, and the central aperture velocity 
dispersion. 

Simulated gas-rich binary galaxy merger remnants show 
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remarkably regular structure when the correct scale con- 
stants are used to plot dimensionless density and kinetic 
energy profiles. The stellar anisotropy is essentially the sole 
source of the variation in the kinematic structure of these 
simulated merger remnants. Dark halo concentration and 
progenitor gas fraction contribute to a systematic difference 
between the Sbc and G series simulations, but within each 
set of simulations the structure is remarkably consistent. 

We emphasize that it is not incorrect to phrase ques- 
tions about homology in terms of r e g. Previous theoretical 
and observational studies that have done so have correctly 
asked and answered well-posed questions about galaxy scal- 
ing relations. However, we find that it is more informative 
to phrase questions about homology in terms of r cq rather 
than r e s, at least for our set of numerical simulations. Using 
r c ff , we would have concluded that our merger remnants are 
not homologous owing to a systematic change in the dark 
matter fraction within r e ff. However, if we instead scale by 
r cq , we find that the merger remnants appear significantly 
more uniform. They are not perfectly homologous, but they 
are closer to homology than we would conclude based on 
r e ff. We are eager to learn whether the same holds true for 
numerical simulations by other groups and whether it holds 
true for actual galaxies. 

In order to facilitate the use of these scaling constants 
to analyse actual galaxies, we calculated the ratios r eq /r e ff 
and Meq/Mbar for a set of observationally plausible two- 
component galaxy models and provided fitting formulas for 
the results. 

The fact that r oq plays such a prominent role in these 
scaling constants immediately recalls the disc-halo conspir- 
acy arising from the study of spiral galaxy rotation curves 



I Bahcall & Casertano 1985 Burstein & Rubin 1985 Kent 



1987). Blumenthal et al. (19861 argued that such a con- 



spiracy will arise naturally if the dissipational collapse of 
baryons is limited by their initial angular momentum and if 
the baryon fraction is roughly equal to the spin parameter 
A = JE 1/2 G~ 1 M~ 5/2 . The result presented here appears to 
be another manifestation of this phenomenon. If the baryons 
have very little angular momentum, they will contract to 
the point where they are completely self-gravitating and the 
dark matter does not play a role in the dynamics, in which 
case r e f; and not r Gq is the important radius. On the other 
hand, if the baryons have too much angular momentum, they 
will not contract enough and the more massive dark matter 
halo will dominate the dynamics while the baryons func- 
tion more as a tracer population. Simulated galaxy merger 
remnants are nearly homologous when r cq rather than r e g 
is used as the dimensional scaling constant, indicating that 
galaxies are in a 'Goldilocks' state where both baryons and 
dark matter play important dynamical roles. This is likely 
because of the initial balance between the baryon fraction 
and the spin parameter of dark matter haloes. 
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